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General-Relativistic Curvature of Pulsar Vortex Structure. 
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Abstract — The motion of a neutron superfluid condensate in a pulsar is 
O ; studied. Several theorems of general-relativistic hydrodynamics are proved for 
^ ' a superfluid. The average density distribution of vortex lines in pulsars and 
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their general-relativistic curvature are derived. 
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1 Introduction. 

m 

O The superfluidity of matter in rotating neutron stars combines with the close- 
3 ness of their gravitational (r^ = 2GM/c^) and geometrical [R) radii, where 
^ {G) is the gravitational constant, (c) is the speed of light, and (M) is the 
neutron-star mass (see, e.g., Manchester and Taylor 1980). Therefore, general- 
; relativistic effects can appreciably affect the processes in the superfluid cores of 
^ ; pulsars and the mechanisms of glitches. In this regard, Andreev et al. (1995) 
>! discussed low angular velocities in general relativity: Q < Qc < Ti/irrfR^)^ 
^ ■ where m* is the mass of a superfluid condensate particle (a Cooper pair), ~ f^c 
^ • (see Kirzhnits and Yudin 1995), and the number of vortex lines (VLs) is and 
1, respectively. 

Here, we deal with the realistic, opposite case, 

Vt,<^Vt< c/R, (1) 

when there is a dense system of VLs. Recently, the interest in this problem 
has risen dramatically, which is most likely attributable to an increase in the 
accuracy of measuring glitches in pulsars. Consequently, it becomes possible 
to detect in principle post-Newton gravimagnetic effects in pulsars, which are 
described below. See, e.g., Prix (2000) and Langlois (2000) for an overview of 

-'^E-mail address for contacts: shatskiy@lukash.asc.rssi.ru , tel. in ASC: (095) 333-3366. 
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this subject. Here, we propose a method of solving the problem in question 
that slightly differs from that proposed by these authors. 

Bekarevich and Khalatnikov (1961) proved that rigid-body rotation is stable 
in the nonrelativistic case, and that VLs are rectilinear and distributed at a 
constant density. Because of general-relativistic gravielectric and gravimagnetic 
effects, which are determined, respectively, by Vgap and rot g (by the definition 
of Landau and Lifshitz (1988b), g = {—goa/ doo} and gik is the metric tensor), 
differential rotation becomes stable (to be more precise, this rotation is stable 
dynamically, while rigid-body rotation is stable kinematically as before; see 
below). In this case, VLs curve and are redistributed in space. These effects 
are described below. 



2 Analyzing changes in the system whenan ordinary fluid 
is replaced with a quantum, fermi superfluid. 

2.1 General Relations for a Superfluid Condensate in General Rel- 
ativity 

The wave function of a superfluid condensate is known (see, e.g., Bekarevich 
and Khalatnikov 1961) to be 

i/j = z/exp(i0). 

Given the identity uiu^ = 10, the generalization of the nonrelativistic relation 
V = to general relativity is 

= ^, e = (2) 

Indeed, let us derive an expression for a general-relativistic superfluid current 
ji = nui f\ with the continuity equation jf = j!- = 0. 

Let us write the Lagrangian of a superfluid condensate in general relativity 
by using the Madelung hydrodynamic representation (see Grib et al. 1980; 
Bogolyubov and Shirkov 1993): 

L = i^*^^^ + l^pFdV^p) = Iz^^^.^c^V + iy'F{u'). (3) 

^The indices i, j, fc, Z, to, n run the series 0, 1, 2, 3; the indices a, f3, 7 run the series 1, 2, 3. 

•^We introduce the notation in which the scalar densities e, uj = p + s, and n are identified with their 
eigenvalues (i.e., with the values in a commoving frame of reference) , respectively: the energy density, pressure, 
thermal function, and density of the particles; Ui are the components of the 4-velocity vector for the matter. 
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The principle of least action relative to a change in (p yields 



= 0. (4) 

Consequently, the current vector is f = z/* = const ■ v'^d'^cj). This is seen from a 
comparison of j* and j*y in the nonrelativistic case: = /c and f = nv'^ jc. 

For a superfluid, we write in this case: n = = is/ ~ m*'c^'^^^'^- 

Since the continuity equation i\ = must hold in any case, we derive Eq. (4) 
for a superfluid. Hence, we obtain for a superfluid current in general relativity 

f = -^ly^d'ct) = nu\ (5) 



or, for covariant quantities, 



Ti 

ji = -—ly^d^cp = nu,. (6) 



m*c 



A scalar multiplication of the latter expression by Eq. (5) yields 

/ n ^2 



Expression (2) follows from this. 

In the Newton approximation, k"^ « {rrfc/hY. The scalar k is identified 
with w/{nch) (see below). Hence, since = n in the Newton approximation, 
Eq. (6) or (2) yields 

h 



m*c 



In seeking to approach an optimum regime, a superfluid current undergoes a 
well-known rearrangement (see Bekarevich and Khalatnikov 1961). As a result, 
while the current remains potential "almost everywhere" , vorticity arises, which 
generally coincides with the vorticity for an ordinary (nonsuperfluid) fluid at 
the same point. This occurs, because a system of VLs is formed, with the phase 
of the wave function = ip + f{t) corresponding to each of them ^. According 
to (7), we then have for a single vortex 

m*c 



Hn what follows, we use a cylindrical coordinate system with the z axis directed along the spin axis: x^'^'^'^ 
ct, p, (j), and z. 
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On the vortex axis itself, where is uncertain, the wave function loses its 
meaning, and the superfluidity vanishes. This corresponds to a physically non- 
superfluid VL "core" of a macroscopically small radius, which is responsible for 
the nonzero curl of velocity (see Bekarevich and Khalatnikov 1961). 

The above reasoning is universal and equally applies to the nonrelativistic 
and general-relativistic cases. 

2.2 The General- Relativistic Bernoulli Theorem 

For the subsequent analysis of a superfluid in general relativity, we need the 
general-relativistic Bernoulli theorem. 

The classical Bernoulli theorem is derived for a steady, isentropic fluid flow 
from the Euler hydrodynamic equations (see Landau and Lifshitz 1988a): 

5V w 
+ (VV)V = -V- 



dt mn 
In the relativistic case, these equations can be written as 

u'di{k*Uj) = djk\ (9) 

where k* = u/n. Note that these equations have a particular solution of the 
form k*Uj = const • dj(f), where is a scalar function of the coordinates. There- 
fore, in contrast to the nonrelativistic case, the quantity k*Uj rather than the 
velocity has the potential. Hence, according to (2), we immediately identify k* 
for a superfluid with chk: 

k* = - = Tick (10) 

n 

In general relativity, dj is replaced with the covariant derivative V^: 

u'V^{k*Uj) = Vjk\ (11) 
or, expanding the covariant derivative, 

u'di{k*Uj) - u'';]k*Um = djk*. 
Substituting the expressions for the Christoffel symbols in this relation, 

Tj = ^P""(5^Pin + djQin " OnQij) , (12) 

yields 



u%{k*Uj) - -k*u'u''djgin = djk*. (13) 
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Denoting 

4)j = dj(f) = kuj] (f)^ = d^(f) = ku^ , (14) 
we obtain for a superfiuid 

(^^5,0, - ^(|)'r^Jg^n = \djk\ (15) 
In that case, since k'^ = gincl)^ (jf- , we derive two equations, 

0"[5«(/)O-5o0a] =0. (16) 

and 

r - d,ct>^] + ctP[d,ct), - d-yM = 0. (17) 
The former foUows from the latter. Denoting 



(7aj = :^[t^a07 - d^(j)a], (18) 



we obtain for cTq,^ f\ 



(^a-f = ^ [0a^70O - (j)jda(j)0 + (j)jdo(j)a - (j)ado(j)j]. (19) 



We emphasize that our system is not axiaUy symmetric because it contains 
vortices, but this symmetry is restored by averaging over the VL distribution. 
As a result, there is no steady state in this case. Clearly, an infinite number of 
stationary frames of references exist for an ordinary fluid uniformly rotating in a 
vessel. Below, we show that there is a (unique) frame of reference for a rotating 
superfiuid comoving with VL cores in which all quantities are stationary. For 
this frame of reference, we derive from (17) dj(f)Q = 15 on all current lines where 
o'a-y = {d^ in any direction). Therefore, 



kuQ = 00 = const. (20) 



'"'Since the condensate ceases to exist on vortex lines, the concept of phase loses its meaning; thus, the 
theorem on the curl of a gradient ceases to be valid on vortex cores, and the 4-gradient of phase is replaced 
with a 4- vector (pi. 

^The prime means that the quantities under consideration refer to a stationary frame of reference outside 
VL cores. 
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This is the general-relativistic analog of the BernoulH theorem. In the nonrel- 
ativistic hmit (see Landau and Lifshitz 1988a): 



1 w m*c. p 
k = T >^—il + ^), Uo 
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hen h pc^ ' \ 1 ~ "^^Z 

^ 1 + + — , =^ A;wo ^ ^ 1 + — ^ = const, 

where and p are the Newton gravitational potential and the matter density, 
respectively. This is the ordinary Bernoulli theorem. 

For a rotating superfluid with vortices, the laboratory frame of reference 
ceases to be stationary, but we can choose a frame of reference for a rotating 
superfluid that comoves with VL cores; in this case, the frame of reference is 
stationary and does not comove with the superfluid. 



3 Determining the relationship between dynamical su- 
perfluid parameters 

3.1 The Principle of Least Action for a Rotating Superfluid in Gen- 
eral Relativity 

Hartl and Sharp (1967) proved that rigid-body rotation with an angular velocity 
equal to the shell angular velocity is most favorable for an ideal, ordinary fluid 
in general relativity. 

Let us prove that this assertion remains also valid for a superfluid Q. 

In view of the general-relativistic Bernoulli theorem and by analogy with the 
study of Hartl and Sharp (1967) for an ordinary fluid, we assume the spatial 
components d^cj) and the probability density u'^ for a condensate particle to be 
detected at a given point to be independent variables for superfluid dynamics. 
In addition, we will remember that Xk are the variables that determine the 
vortex shape and coordinates, while the metric components gik are the variables 
that determine the system's gravitational field. 

Given the conservation of momentum and the total number of particles, the 
expression for the total energy of a superfluid can be written as (Hartl and 

^By the angular velocity of a superfluid, we imply the mean value of u'^ /{cvP) averaged over the surface 
orthogonal to the vortex direction near the point under consideration whose area is much smaller than the 
system's cross-sectional area, but, at the same time, is still much larger than the square of the mean separation 
between vortices. As we show below, this condition can definitely be satisfied for pulsars. 
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Sharp 1967) 



1 = zy20Q0O _ ^,2 1^ ^ ^2/2] - zy^Vcp) - /^^'/, (21) 

c 

where Q and /j, are the corresponding Lagrange factors. 

In addition, we must take into account the Lagrangian Lq of the shell. Given 
momentum conservation, the analog of (21) for the shell can be written as 

where Jo and Qq are the moment of inertia and angular velocity of the shell in 
the frame of reference under consideration, respectively. Varying over r^o leads 
to the equality 

n = Qo- 

If, however, the satisfaction of the general-relativistic potentiality condition for 
a superfiuid flow must also be taken into account, then, according to the general 
rules for imposing additional conditions, we must add the term AL to L to 
express the fact that the circulation of the phase gradient over a closed contour 
is proportional to the number of vortices crossing this contour. This is the 
general-relativistic generalization of the potentiality condition for a superfiuid 
flow: 

f d^cf) dx^ = j {dpd^cp - d^dpcp) dx^ A dx^ = 27rK, (22) 

where dx^ A dx'^ is the directed surface element pulled over contour F, and K 
is the number of vortices crossing the contour. Thus, the general-relativistic 
potentiality condition, according to (18), can be written as 

N 

27ra/3^ = E 27r4^^'(x - xu), (23) 

k=l 

while the addition to L corresponding to this condition can be written as 

N 

AL = ^('\x)[27rap^ - E ^i^n^.S^x - Xk)]. (24) 

k=l 

Here, according to (18), 27r(7/3^ is the curl of 0^; n^^ is a unit antisymmetric 
tensor dual to one of the surfaces, which are orthogonal to the vectors of the 
vortex direction at each point, n^^ is defined at the point of intersection of 
this surface with vortex /c, and the summation over k is performed over all 
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vortices in the system 5'^{x — Xk) is the bivariate delta function, and x is the 
radius vector defined on this surface; and the antisymmetric tensor ^^^{x) is 
the Lagrange functional factor. 

4 Determining the Relationship between the Compo- 
nents of the 4-Gradient of Phase for a Superfluid Con- 
densate 

Action A in general relativity is known to be related to Lagrangian L by 

A = J LdX, dX = y/^d^X, 

where the integration is performed over 4-space; g is the determinant of the 
metric tensor gik. We denote the integral of L over space by 

E = J ^/^Ld^x, i^^ = J V^^Ld^x). (25) 

By its physical meaning, E is the system's energy for imposed additional con- 
ditions, such as allowance for the conservation of angular momentum, the total 
number of particles in the system, etc. 

We break up the integral over space into an integral over the surface dual 
to the tensor n^^ and integrals over the lengths of the vortices orthogonally 
crossing this surface. Because of the presence of delta functions, the following 
sum remains from AE: 

r ^ 

^E= ^g2e\x){dp(t),)d'x + E 2i,hV^ge\xk)nl^, 

k=i 

where Ik is the length of vortex line k. 

After varying AE over 0^, the following term remains: 

- 2dp{^gei (26) 

If we discard addition (26), then the subsequent analysis will be valid only 
for an averaged description of the superfluid motion. Let us prove that term 
(26) vanishes on vortex cores. Since the coordinates Xk must correspond to 
equilibrium vortex positions in the system, the system must be stable against 
core displacements orthogonal to the direction of the vortices themselves: 

'^^kn^. m^gem^^., = o. (2?) 

^In general, addition (24) to L should be written for each such surface, but this is implied by default. 
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This proves the above assertion. The corrections related to the term in ac- 
tion will not be considered everywhere, because we are interested in the system's 
dynamics only when expression (27) holds, i.e., on vortex cores, or, equivalently, 
only an averaged description of all quantities for the system. Below, we thus 
denote average quantities by a hat above the symbol (e.g.. A). According to 
Eq. (3), the principle of least action relative to a change in i/'^ leads to the 
equation 

I4 = #2{-'[W2 + F(.^)]} = 0. (28) 
Given that k'^ = g^^(j)a4>i3 + '^9^^4>Q4>a + 9^^{4>qY and denoting 

A" = B- = If, (29) 

we obtain 

= 2p"^0/3 + 2pO"0o + 2p°"0aA" + 2pOVo^" = 2[0" + /A"]. (30) 



Taking a variational derivative of Eq. (21) with respect to (f)^ and z^^, using 
Eqs. (28), (29), (30), we derive for the average quantities 

B^i4>o + - /^) = 0^ - ^4>X- (31) 

^0 + -0^ - /i = 0. (32) 
c 

Comparing Eqs. (31) and (32), we obtain the analog of rigid-body rotation for 
a superfluid: 




4.1 Magnus Force and the General- Relativistic Theorem on the 
Conservation of Circulation 

The Magnus force acts on a rotating body in an incoming flow and is at- 
tributable to a nonzero pressure difference for the opposite sides of the flow 
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around the body. In tern, the pressure along the body boundary changes be- 
cause of the Bernoulh theorem: the velocity of the medium that flows around 
a rotating body changes when going around the body axis. 

Let us write the Bernoulli equation for the nonrelativistic case: 

p + = const (34) 

It would be natural to consider a portion of the VL core as the body on which 
the Magnus force acts. As we show below, the Magnus force does not depend 
on the radius a of this core. Choose a cylindrical coordinate system whose z 
axis coincides with the rotation axis of this core. Since the core radius a is much 
smaller than any scales on which the incoming flow produced by the remaining 
vortices changes appreciably, this flow may be considered constant for the flow 
around the core. According to (33), the physical velocity of this flow at the 
core point is 

Vo = [fi X r] 

At the same time, the velocity produced by the core itself on its boundary, 
according to (8) is 

Vfc = ^^[^z X a], 

where a and are the vector in the direction of the core radius (equal to a 
in magnitude) and a unit vector along the core rotation axis, respectively. The 
total velocity on the core boundary is 

^ = ^o + ^^[e. xa]. 

The force per unit core area is equal to the pressure multiplied by a unit vector: 
— Ba = —a/a. Accordingly, the force per unit core length is 

27r 

¥ = - j eap{(p) dl = - j eapi(p)a dip. (35) 



Expressing p from Eq. (34) and substituting it in Eq. (35) yield 

27r 

F = - f a{const - mn(Vg + + 2VoVfc) /2}d<f. (36) 



Since the (p — independent terms vanish and since n and m ^ m* for a 
superfluid in the nonrelativistic case, we derive for the Magnus force per unit 
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core length 




It follows from Eq. (37) that the Magnus force acting on the VL core that is at 
rest relative to a remote observer causes it to move toward the vessel wall. As 
it accelerates toward the wall, an incoming flow emerges (to be more precise, 
the core itself runs on the superfluid); as a result, the Magnus force changes its 
direction, causing the VL core to precess around some point of the superfluid 
flow. This precession is rapidly damped, and the core starts moving in such 
a way that the Magnus force does not act on it, i.e., that the system's energy 
becomes minimal. This requires that the core be at rest with respect to the 
superfluid flow at its location. The above reasoning proves that the VL system 
is frozen in, i.e., that there is no slip in the nonrelativistic case. 

In general relativity, expression (34) for the Bernoulli theorem is replaced 
with expression (20). 

For the nonrelativistic case, Hess (1967) showed that the VLs are, as it were, 
frozen in a superfluid — move at angular velocity Q, the shell rotation velocity. 
Thus, there is no slip relative to this angular velocity. 

To generalize the slip theorem to general relativity, we do not need to repeat 
similar calculations in order to determine the Magnus force, which, incidentally, 
are very complex. It will suffice to note that 0o = kuQ, k = /c(|V|, j9, n, ^^ifc), 
and uq = uo{\V\, gik) and that when going around the core of a vortex line, the 
changes in metric gik and density n, if any, are so negligibly small f\ that they 
may be disregarded. 

Therefore, we can write for the superfluid portion adjacent to the core 

p = p{\V\) (38) 

We thus see that there is no general-relativistic Magnus force for a VL motion 
with |V| = const around the core. This is possible only when the VL core 
accompanies the superfluid flow, i.e., when V is produced by the core itself. 

Consequently, given Eq. (33), we see that there is no slip in general relativity 
either. 

In conclusion, note that the absence of slip also follows from another impor- 
tant theorem of hydrodynamics, the theorem on the conservation of circulation 
(see, e.g.. Landau and Lifshitz 1988a). 

^The core radius is known to have sizes of the order of 1 10 interparticle separations. 
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5 Calculating the mean density and curvature of vortex 
lines in a pulsar with general-relativist ic corrections 

5.1 Passing to a Rotating Frame of Reference 

As will be evident below, it is more convenient to perform an analysis in a 
comoving (with vortex cores), i.e., rotating frame of reference. 

According to Landau and Lifshitz (1988b), the following formulas can be 
derived that relate the tensor components in various frames of reference: 

9pp ~ 9pp'i Qifip ~ 9<ff'i 9zz ~ 9zz 
' 9o^ = 9ov + ^9w; pSo = Poo + (?)^Pw + 2|po^ (39) 
u'p = Up-, u'^ = u^; u'^ = Uz] Uq = uo + -u^. 

Here, the components in a frame of references rotating with angular velocity Q 
are marked by primes. 



5.2 Calculating the Covariant Curl of Superfluid Velocity and the 
Vortex Density in the System 

By definition, the mean density of vortices is the number of vortices crossing 
the orthogonal surface divided by its area. Therefore, using Eq. (22) to derive 
the vortex density crp^, we obtain 

It is convenient to express the quantities (J)^ in terms of 0o, the metric, and the 
shell angular velocity in the frame of reference under consideration: 

(pa = (P^90a + (P^9aj = /(POa + 9aj^), 
00 = /poo + </>^P07 = (p^i9oo + 90j^)- 

Hence, according to (33), we obtain for the average quantities 

(pa = 00 • (41) 

9oo + -9(pO 

Since the average quantities do not depend on time and angle ip, we introduce 
the notation 

= 27r<7^^ = -5^0^; y = 0o; z = (42) 

fi'oo + -^9ipO 
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and derive from (41) 

X^ = d^iYZ). (43) 
On the other hand, we obtain from Eqs. (19) and (33) 

= {c/Q)d^Y, (44) 

Solving the last two equations for and Y yields 

X^ = const ^-^^Jn^y . (45) 



c 

const , 

y = (46) 



Since Y = ^ k ^ nfc/Ti in the nonrelativistic limit, we see that const = 
nfc/Ti. Hence, we have 

As was shown above, il' = in the comoving (with cores) frame of reference; 
therefore, Eq. (47) in this frame of references is especially simple: 



a' _ .a.^ 



d,9' (48) 



where, according to Landau and Lifshitz (1988b), Qr^ = —Qo-y/goo, 7 = —g/goo is 
the determinant of the spatial metric tensor, and the three-dimensional vector 
a^lpha dual to the tensor cr^/j was defined as cr" = (2^^)"^ • e'^^^a^p, c^^^ is a 
unit antisymmetric tensor. The vector coincides in direction with the vortex 
direction in the system and is equal in magnitude to the mean vortex density 
at a given point. 

For the invariant mean vortex density, we can write 



a = ^a.ja^J = sj2a^a^^- + dapa^^ (49) 

As a result, we have for the invariant density in the first post-Newton approxi- 
mation 

^^\^^a\y/¥^. (50) 
which matches the nonrelativistic limit for the VL density. 
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According to Eqs. (39), we have for g'^ 



9^ 



a 



goo -"fi 



^ C. ^ nnr\ ^ 

Given that g' = g, we obtain for 7': 



-f 9' 7 



a" ^ = (5?^^ = const. 



^ 900 l + ©^|f-2?^o.- 

In the nonrelativistic case, 7 p^, goo 1, g^^p — p^, gzz —1, go-y 0, 
and 4)0 ^ k ^ m*c/h; therefore, we derive the already known expression in 
this hmit 

For a milhsecond pulsar, (Jq ~ 10^^ cm^^, with a separation between vortices 
d ~ 10^^ cm corresponding to this value. As for the radius of the core itself, its 
order-of-magnitude value is 10~^^ cm. 

As we see from this section, the general-relativistic corrections that deter- 
mine the curvature of VLs and the change in their mean density are small 
for real pulsars, and, therefore, the relative curvature does not exceed a few 
percent. 



6 Calculating corrections for a homogeneous model. 

Let us calculate the invariant density of vortex lines in a pulsar in the first 
post-Newton approximation. The model is based on the assumption that the 
pulsar interiors rotate at angular velocity and, because the compressibility of 
a neutron condensate is low, its density is assumed to be p in the entire volume 
0. 

In the first post-Newton approximation, the metric in a conformally Eu- 
clidean coordinate system Q can be written as 

ds'^ = (1 + 2(p)dt^ - (1 - 2(^)[(i/?^ + dz^ + p'^dip'^] - 2g^dtdip. (52) 

According to Eq. (50), we write 

a = Vr^^^7 K - ^9^ J 900] , ('^0 = ^) • (^3) 



^°For convenience of calculations, we take G — 1 and c = 1 in this section.) 

^^It is easy to see that the result does not depend on the choice of a coordinate system in this approximation.) 
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Since g"^ = 1/gaa-, the corrections to gq can be arbitrarily divided up into two 
groups: 

(1) Gravimagnetic corrections: 

m* {dp + d^) 

= w — -p — • 

(2) Gravielectric corrections: 

We thus see that the VL curvature and redistribution in a pulsar result from 
the gravimagnetic interaction of VLs with the Lense-Tirring field of the system 
and from the gravielectric deformation of the system's Euclidean geometry. 
The Newton gravitational potential (p of the model can be easily calculated: 

(p = -27TpR\l - x^S - y^S) (x = p/R, y = z/R). (54) 

Hence, it is easy to derive an expression for the gravielectric corrections: 

(72 = 47rpi?^(7o(l -x^ - y^jZ - 2xy/2>). 

Given that (p^^ = —AirpR'^/S on the stellar surface, we obtain for the gravielectric 
corrections 

a2{x, y) = 3|(^^|(7o(l - x^ - y^ jZ - 2xy/'i). (55) 

To determine the gravimagnetic corrections, we use Eq. (106.15) from Landau 
and Lifshitz (1988b) to derive the metric components g^. In Cartesian coordi- 
nates. 



V 



where Ua = (r^ - r^)/|r - r'|. 

Hence, it is easy to calculate the angular component of the gravimagnetic 
field in cylindrical coordinates: 



/ N P f ~ 7 i+TT 7 r / 7 ,r7^p'cosp r^pV cos a cos 7 , 

9^^) = -\ J pdp\tldz'p'dp'{ + ^1 _ (57) 

z y |r r I |r r I 

Here, (r — r')^ = {z — z')^ + J>^, = p^ + p>^ — 2ppl cos p, 

cos a = cos(^sin^ + sin(y5C0S^, cos 7 = — sin^, sin^ = ^sin(/?, and 
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2,2 /2 

COS e = p Y"^ • 

As a result, expression (57) reduces to 

/■ , 7p'cosip p' sinV(p'^ + PP - PP' cos V?) . 

^v'(r) = y p^i^^ll^c^^ p dp { + 1 H3 }• 

^ y I* * I I* * I 

(58) 

Since the integrand is even in variable by changing variables: x = p/R^ 
x' = p' / R, y = z / R, y'{y) = {z' — z)/R, and p'{x, if) = p/R, we can write the 
expression for the gravimagnetic corrections as 



/I J- 

ai{x, y) = ao\(fJ I dip I dx'f{x, x', y, (p) (59) 



where 

f(x, x\ y, ip) = -^Ad:, + dy) / dy' (r^x^ + ^^^^^"'^(f;+p;^7f^'^°^^) | , 



V'\ = -Vl - a:^'^ - y, y'2 = +Vl -x''^ -y. 

(60) 

Hence, integrating and then differentiating yields 
A = 7x cos 

5 = icx' sin^ if{x''^ + p'a: — p'ar' cos if) , 
C = Vy'" + p'" 

= 7 COS 

= sin^ ip{x''^ + p'ar — p'a:' cos if) + xar' sin^ ifijp' -\- {x — x' cos iff' j pi)-, 
4 = A, ln(^' + C) + ^^g;f|f + ^ - By'^^^ - 25y'^=^, 

/y = -A/C-B/C3, 

^!{x,x\yM = -ti\^.^ht\ 



When deriving the last expression, we took into account the fact that the 
derivative with respect to y could be taken inside the integral and that dy = 

-dy.. 

Integral (59) can be calculated numerically. The integration results with 
gravielectric corrections (55) are shown in the figure (the maximum amplitude 
of relative corrections is 1). 



16 



Fig. 1: Dependence of corrections (T2) I {^o\'f> ^\) along the vertical axis on coordinates: 
the rotation axis and the perpendicular axis in the equatorial plane (the lower right 
quadrant of the pulsar meridional section is shown). 
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